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We revisit classical nucleation theory (CNT) for the homogeneous bubble nucleation rate and 
improve the classical formula using a new prefactor in the nucleation rate. Most of the previous 
theoretical studies have used the constant prefactor determined by the bubble growth due to the 
evaporation process from the bubble surface. However, the growth of bubbles is also regulated by the 
thermal conduction, the viscosity, and the inertia of liquid motion. These effects can decrease the 
prefactor significantly, especially when the liquid pressure is much smaller than the equilibrium one. 
The deviation in the nucleation rate between the improved formula and the CNT can be as large 
as several orders of magnitude. Our improved, accurate prefactor and recent advances in molecular 
dynamics simulations and laboratory experiments for argon bubble nucleation enable us to precisely 
constrain the free energy barrier for bubble nucleation. Assuming the correction to the CNT free 
energy is of the functional form suggested by Tolman, the precise evaluations of the free energy 
barriers suggest the Tolman length is ~ O.Scr independently of the temperature for argon bubble 
nucleation, where a is the unit length of the Lenard-Jones potential. With this Tolman correction 
and our new prefactor one gets accurate bubble nucleation rate predictions in the parameter range 
probed by current experiments and molecular dynamics simulations. 


I. INTRODUCTION 

Bubble nucleation in liquid is a liquid-to-vapor tran¬ 
sition [l[ phenomenon, and plays an important role in 
many areas of science and technology, e.g., degassifica- 
tion of steel [I|, bubble jet printers vulcanism 
the direct detection of dark matter [6l-l9| and medicine 
©[HI- Bubble nucleation occurs in a metastable liquid 
under a pressure below its equilibrium vapor pressure. 
The liquid phase can even exist at negative pressures, and 
bubble nucleation is often observed at negative pressure, 
although there are no fundamental differences between 
a boiling, superheated (positive pressure) and a cavitat- 
ing, stretched (negative pressure) liquid [l^- Studies of 
homogeneous liquid-vapour nucleation typically use the 
classical nucleation theory (CNT) for the bubble nucle¬ 
ation rate. However, the range of applicability of the 
CNT is not well understood. 

Numerical techniques such as molecular dynamics and 
Monte-Carlo simulations are powerful methods which can 
resolve details of the nucleation process and provide use¬ 
ful test cases for nucleation models [isl - l^ . Typically, 
these simulations show large deviations from the CNT 
predictions. The CNT nucleation rates usually underes¬ 
timate the bubble nucleation rates by very large factors 
[Tsl - fisl . [20l - [2^ m, [ 2 ^ . Most of the simulations for bub¬ 
ble nucleation in the literature use around 10^ or fewer 
atoms, making it difficult to measure nucleation rates di¬ 
rectly. 

Recently, Diemand et al. [ 2 ^ presented large-scale, 
micro-canonical molecular dynamics simulations of ho¬ 
mogeneous bubble nucleation with 5 x 10® Lennard-Jones 
atoms, and succeeded in measuring nucleation rates in 
the range of 10^^“^®cm“®s“^ for argon by directly resolv¬ 


ing bubble nucleation events in the steady state nucle¬ 
ation phase. The measured rates agree well with the CNT 
within two orders of magnitude in the superheated boil¬ 
ing regime (positive ambient pressure), while the CNT 
prediction underestimates the nucleation rates signifi¬ 
cantly in the cavitation regime (lower temperatures and 
negative pressures). 

The kinetics of explosive cavitation in liquid has also 
been investigated in laboratory experiments, which have 
measured the superheat temperature for liquid argon at 
both positive and negative pressures, by pulse hea ting 
liquid around a thin wire in a negative pressure wave |l2l . 
[ 2 ^. Using this method, Vinogradov et al. [l^ measured 
nucleation rates of 10^®“^®cm“®s“^ in superheated liquid 
argon of high purity. 

The recent advances in molecular dynamics simula¬ 
tions and experiments enable us to precisely test the¬ 
oretical bubble nucleation models and also to improve 
them. One of the most serious problems in the CNT 
model is that the bulk value of the surface tension is used 
to evaluate nanobubble formation energy. Since the nu¬ 
cleation rate of the CNT depends exponentially on the 
formation energy, an incorrect estimate of it can cause 
a huge error in the nucleation rate. According to Tel¬ 
man’s correction [s^, the surface tension at surfaces of 
small nuclei (bubbles or droplets) is dependent on their 
radius. A model parameter called the Tolman length can 
be determined by the measurement of the surface tension 
of small nuclei [3ll - l40l| . Diemand et al. [ 2 ^ showed 
that the CNT model agrees well with the nucleation rate 
from their MD simulations, by using such a model for 
the surface tension with a proper Tolman length. The 
introduction of a Tolman correction can also significantly 
improve model predictions in vapor-to-liquid droplet nu- 
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cleation [40l |. 

In bubble nucleation, additional detail in the treat¬ 
ment of the process is required, in comparison to droplet 
nucleation, because of the variable vapor pressure and 
density in bubbles as they grow . The vapor pres¬ 

sure in bubbles varies as they grow and has a significant 
effect on their growth rates and on the pre-exponential 
factor in the CNT formula for the nucleation rate. Thus 
we have to solve the two-dimensional evolution (i.e., the 
radius and pressure) for the bubble growth process. An 
extensive study of this problem has already been done by 
Kagan [4l| . Kagan showed that the pre-exponential fac¬ 
tor is strongly dependent on the ambient liquid pressure 
reduced from the saturation. Despite this, a constant 
prefactor is usually adopted in the widely-used CNT. Fur¬ 
thermore, the formation energy of a bubble also depends 
on the vapor pressure. Although the vapor pressure is 
approximately given by the saturated pressure for criti¬ 
cal bubbles, a more accurate vapor pressure is necessa^ 
for the evaluation of the formation energy in the CNT [J] . 

Although these treatments have been developed in in¬ 
dividual studies, they are not included directly in the 
widely-used CNT. For comparisons with recent molecu¬ 
lar simulations and laboratory experiments, we should 
use a precisely crafted expression for the CNT bubble 
nucleation rate. Detailed comparisons also enable us to 
correctly determine the Tolman length in the model of 
the surface tension. 

In this paper, we first present a more complete expres¬ 
sion for the nucleation rate, by summarizing the above 
studies (Sec. II). Next, we compare the improved model 
with the original CNT or with the measured values in 
the recent MD simulations and experiments (Sec. III). 
We find that the difference between the improved model 
and the CNT could be several orders of magnitude for re¬ 
alistic bubble nucleation parameters. Comparisons with 
MD simulations can determine the Tolman length more 
accurately thanks to the more accurate pre-factor in our 
model. From our comparisons at various temperatures, 
the Tolman length is obtained as ~ 0.3 cr, where cr is the 
unit length of the Lenard-Jones potential. A summary 
of our findings can be found in Section IV. 


II. NUCLEATION THEORY 


where T is the temperature, k is the Boltzmann constant, 
AG{i) is the minimum work for the formation of a bubble 
with i molecules and no is the number density of liquid 
molecules. 

In the classical nucleation theory (CNT), the minimum 
work for the formation of the critical bubble is 


AGcnt 


16717^ 

3(Peq-Pl)2’ 


( 3 ) 


where 7 is the surface tension, Peq and Pi are the equilib¬ 
rium vapor pressure at saturation and the liquid pressure, 
respectively. The prefactor Jq is given in the CNT by 


do,CNT = 



( 4 ) 


which is the same as in the droplet nucleation. In the 
above, m is the molecular mass. Thus, the widely-used 
expression of the CNT nucleation rate is 
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We will describe more accurate expressions for AG and 
Jo in subsections B and C. 


B. Free energy for bubble formation and the 
Poynting correction 

The minimum work AG for the formation of a bubble 
with radius r is given by [l|, 

AG [Mg(^g) - Mm - {Pe - fl)} 

-p 47rr^7, (6) 

fig and fii are the chemical potentials of the gas and liquid 
respectively, Pg is the gaseous pressure in the bubble, and 
the molar volume of a gas Vg is given by kT /Pg, assuming 
the ideal gas. Eq. dH]) also assumes spherical bubbles. 
Note that the work for bubble formation is a function of 
two variables (r, Pg) or (r, z). The number of molecules 
in the bubble is given by 


A. Classical expression 


Anr^Pg 

3kT 


( 7 ) 


The bubble nucleation rate is the number of stable 
bubbles formed per unit time per unit volume and is 
given by P, 

J = Jone(Zc), (1) 


where ic is the number of gaseous molecule in a critical 
bubble, Jo is the prefactor in the nucleation rate, and the 
number density of bubbles ne{i) is given by 


ne(z) = no exp 


AG(z) \ 
kT J ’ 


( 2 ) 


Previous studies have investigated the minimum work 
for a bubble formation with two variables, taking into 
account bubble compressibility (dll - ld^ . Here we use the 
two variables r and Pg. 

The work AG has a maximum value along the path of 
the bubble growth. Such a maximum point on the path 
corresponds to the critical bubble and is given by the 
saddle point in the two dimensional plane (r, Pg) [II]. In 
the next subsection, we describe the growing path in the 
vicinity of the saddle point. The radius and the internal 
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pressure of the critical bubble are thus obtained from the 
conditions: 


dPJr 


/dAG\ 


= 0 . 


( 8 ) 


Since the first-order derivatives of the minimum work are 
determined by 


5AG\ 

)p 


47rr^ 


kT 


[Mg(^g) Ml] 
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and 


dAG\ 

dPJr 


47rr^ 
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the critical size of the bubble and the critical gas pres¬ 
sure Pg_c are determined by 


Mg('fg,c) — Mb 
27 


( 11 ) 

( 12 ) 


Eqs. (HU and m indicate the chemical equilibrium and 
the mechanical equilibrium at the saddle point, respec¬ 
tively. Using Eqs. (HU and HU, the minimum work for 
the formation of the critical bubble is 


AG (re, Pg,c) — 




(13) 


Using Eqs. (|U-(HU, the second-order derivatives of the 
minimum work at the saddle point are given by 


/ d^AG \ 

\ ) F 

/d^AG\ 

V dPi A 


= —Sttj < 0 , 


3^ 


> 0 , 


and 


/d‘^AG\ 


g,C 


= 0 , 


(14) 

(15) 

(16) 


With Eqs. (IT^ - lfTO . the work for bubble formation 
around the critical size is given by 

AG = - 4.,(r - r^r + P.,c)^(I7) 


up to second order accuracy in (r — Tc) and (Pg — Pg,c). 
This expression shows that the point (rc,Pg,c) is indeed 
the saddle point. 

The gas pressure, Pg,c, in the critical bubble is de¬ 
termined by Equation HU- Integrating = vdP, we 
obtain 


Mg(-fg.c) Mg(-feq) — kTln{ ), 

-^eq 


for gas and 


Ml(-Pl) ~ Ml(-Peq) = ^ll(Pl — Peq), (19) 

for liquid. In Eq. (IT51) we used the equation of state for 
ideal gas, Vg = kT/Pg and vi is assumed to be constant 
in Eq. ((^. Noting fJ.g{Peq) = M/(-Peq), Eqs. (ITT|), (HH) 
and (fTOll yields 


In 


p 

^g,C 

Peq 



( 20 ) 


In (1^01) . Peq(= kT/Peq) is the molecular volume of ideal 
gas at the equilibrium pressure. 

In the classical theory, the right hand side of Eq. (EU 
is set to be zero because the ratio vi/vg is small. This 
approximation in Eq. (17) gives Pg_c = Peq- Then, in the 
CNT, the critical radius is given by 


r-c.CNT = _ - ( 21 ) 

-^eq -^1 

Blander and Katz [l| keep the small term proportional 
to v\/veq in Eq. (1201) and derive more accurate expres¬ 
sions for the critical radius and the nucleation rate than 
CNT. They termed this correction to CNT as the Poynt- 
ing correction. Here we label this correction as PCNT. 
In PCNT, the gaseous pressure in the bubble is obtained 
as 


= exp f 

Peq \ "^^eq Peq / 

~ 1 AL i ^1 -Pi - -Peq 

^eq \ Peq ) 2 \ 'Ceq Peq 


, ( 22 ) 


then we have 


Pgp - Pi = (Peq - Pl)^, (23) 

where 5 is the Poynting correction factor given by 


II 

1-AL 


f—)" 


Veq 

2Peq 

\VeqJ 


The expansion in Eq. (1221) is valid even if P\ is far from 
Peq because of the factor v\/veq is small (^ 0.01 — 0.1). 
This indicates that Pg{rc) is close to Peq even if Pi is 
far from Peq. But, the gaseous pressure Pg can deviate 
considerably from Peq(r) when r is far from re. We find 
that in the Poynting correction factor derived by Blander 
and Katz, the second order term of the small ratio vi/veq 
is incorrect. 

Note that the approximation of ideal gas is inaccu¬ 
rate at a high temperature where the equilibrium pres¬ 
sure Peq(T) is considerably large. We also derive the the 
Poynting correction factor for such a non-ideal case with 
the first order accuracy. Since Pg_c is close to Peq for 
critical bubbles, Vg^s can be replaced by feq(= 'Cgas(Peq)) 
in the first order approximation and we obtain instead of 
Eq.dUl) 


( 18 ) 


Mg(-fg,c) Mg(^eq) — 'Ceq(Pg,c ^eq)- (25) 
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As for Weq, we include the non-ideal effect. Using the 
second virial coefficient B^. the molecular volume Ueq of 
non-ideal gas is given bv|^ 


v,^ = kT/P,^ + B2{T). (26) 


Although we use Eq. (l25l) instead of (18), we obtain the 
same Pg_c as Eq. (22) up to the first order term and thus 
have 5 = 1 — v\/veq- Hence, by using Veq of Eq. (l26ll in 
Eq.(23), we can obtain the Poynting correction factor for 
non-ideal gas with first order accuracy. 

The critical radius in the PCNT is obtained as 


?'c,PCNT = 


27 

(Peq-Pl)S' 


(27) 


With this Poynting correction for Tc, the nucleation rate 
is given by 


J = Jong exp 


IffTry^ 

3kT(P,q - Pi)252 


(28) 


We will consider the prefactor Jq in the next subsection. 
Although 1 — 5 (~ m/ug) is usually small, the difference 
in the nucleation rate between the CNT and PCNT can 
be large because of the strong exponential dependence. 
This correction is necessary especially at a relatively high 
temperature where the ratio vi/vg is not so small due 
to high Peq- At such a high saturated pressure, it is 
also necessary to include the non-ideal gas effect for Veq 
(Eq. [26]). 

For a more accurate evaluation of the nucleation rate, 
we also need to take into account the deviation of the sur¬ 
face tension from the bulk value for the nano-sized critical 
bubbles. In Sec. Ill, we adopt the Tolman correction for 
the surface tension [27] and fix the model parameter, the 
Tolman length, by using the results of MD simulations. 

Before proceeding to the evaluation of the prefactor, 
we show some examples of AG in Figure |T| For the com¬ 
parison with molecular dynamics simulations in Section 
3, here we consider a Lennard-Jones liquid. In Figure |TJ 
panels a and b show AG as a function of the bubble 
radius r for the temperatures T*(= kT/e) = 0.855 and 
0.7, respectively, where e is the binding energy of the 
Lennard-Jones potential. At each temperature, the liq¬ 
uid pressures are set to be 0.017 ecr“^ and -0.16 ecr“^, 
respectively, where a is the unit length in the Lennard- 
Jones potential. The equilibrium pressures are given by 
0.046£cr“^ and O.OlOecr”^ and the surface tensions are 
0.089 £o-“^ and 0.33 £(t“^ respectively at each tempera¬ 
ture [1^. Here, we set vi/vg = 0. To evaluate AG, 
we also fix i (or Pg). In Panels a and b, we plot AG 
for bubbles in mechanical equilibrium (solid lines) and in 
chemical equilibrium (dashed lines). The thin solid lines 
show AG for various constant bubble molecule numbers 
i. Panel c and d show the relation of i and r for the 
mechanical and the chemical equilibriums at T* = 0.855 
and 0.7. 

In the positive pressure case of Panel a, both AG in 
the mechanical and the chemical equilibria are similar 
to one another. They agree exactly at their maxima - 



FIG. 1. (Color online) Free energy for bubble formation AG 
as a function of the bubble radius for Lennard-Jones system 
with T* = 0.855 (panel a) and 0.7 (b). The liquid pressures 
are Pi — 0.017£cr“^ and Pi = —0.16£cr~® for T* = 0.855 and 
0.7, respectively (see text for detail). The black solid curves 
assume the mechanical equilibrium, while the blue dashed in¬ 
dicate the chemical equilibrium. The thin solid curves show 
AG for various constant bubble molecule numbers i. Panel c 
and d show the relations between i and r in the mechanical 
and the chemical equilibriums at T* = 0.855 (panel c) and 
0.7 (d). In the positive pressure case with T* = 0.855, the 
paths and AGs in the two equilibriums are close each other, 
whereas the deviations in them are large in the negative pres¬ 
sure case. Panel d also shows the path of bubble growth across 
the critical size by the doted line. 


corresponding to the critical size. At other radii, AG is 
slightly smaller at chemical equilibrium. The minima of 
each constant-f line are located on the line of the me¬ 
chanical equilibrium. Thus we find that the critical size 
corresponds to the saddle point. 

In the negative pressure case of Panels b and d, the 
mechanical equilibrium lines deviate considerably from 
those in chemical equilibrium. The maximum of AG in 
the chemical equilibrium corresponds to the minimum of 
the mechanical equilibrium case. This also shows that 
it is the saddle point. Panel d shows that the num¬ 
ber of molecules of growing critical bubbles increases for 
the chemical equilibrium case whereas it decreases for 
the mechanical equilibrium case. These facts indicate 
that chemical equilibrium is more realistic than mechan¬ 
ical equilibrium. In the next subsection, we examine the 
growth of bubbles using the model of Kagan [4l| to ob¬ 
tain the prefactor. It also gives us the path around the 
critical size. In Panel d, we also plot the obtained path 
across the critical size (dotted line), which is close to the 
line of the chemical equilibrium. 
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C. The nucleation rate prefactor 


Compared to droplet nucleation, bubble nucleation re¬ 
quires a more detailed treatment due to the bubble com¬ 
pressibility. The vapor pressure in the bubbles varies as 
they grow and significantly affects their growth rates and 
the pre-exponential factor in the CNT expression for the 
nucleation rate, as well as the bubble formation energy. 
Kagan solved the two-dimensional bubble evolution 
and showed that the pre-exponential factor is strongly 
dependent on the ambient liquid pressure. Based on the 
method of Kagan [dlj, we evaluate the exact prefactor 
Jq. 

The nucleation rate is usually given by 

^ (29) 

where the Zeldovich factor Z is 


Z = 


1 fcPAG\ 
2TrkT I di"^ ) ■ 

\ ' Ic 


(30) 


and D is the diffusion coefficient in the i-space given 
by mi 


where a is the evaporation coefficient (often taken to be 
unity). The factor e is introduced to include the effect 
of the temperature difference between the ambient liquid 
temperature T and that at the bubble surface [dll which 
is given by e = adqrc/{XV2'KmkT), where q is the la¬ 
tent heat of evaporation per molecule, A is the thermal 
conductivity coefficient of the liquid, and d = dPgq/dT. 

Bubbles smaller than the critical size shrink, while 
larger ones grow. At the critical size, the growth rate 
therefore vanishes. So we can put dr/dt around r ~ rc 
as 


J = A(.-..), (37) 

Since the vapor pressure in the critical bubble is approx¬ 
imately equal to Peq, the vapor pressure is given by 

dP 

Pg = Peq + ^(r-r,), (38) 

around r ~ Tc, where we have neglected the small terms 
proportional to vi/vg. 

Substituting Eq. (l38l) into Eq. (l36l) and transforming 
the left hand side of Eq. ((551) by the use of Eqs.([7]), (1571) . 
and (1551) . we obtain 


and Dc = D{ic). Note that Z and D^. are evaluated at 
the critical size ic (or Tc). Since the denominator and 
fraction in Eg. (1311) both vanish at the critical size, we 
evaluate it with the second derivatives 


dr 


1 + e 

aVth 


PeqA 



1 + e 
3aEth 



-1 


(39) 


where Vth is the thermal velocity defined by \/kT /(27rm). 
Substituting Eqs. (l37l) - (159l) into Eq. (1551) . we obtain the 
equation for A: 


"—Id 


d^AG 

dP 


(32) 


From Eq. (l29l) . the prefactor in the nucleation rate is 
Jo = ZDc- We use r instead of i in the equations for 
convenience, and rewrite Jq as 



d^AG 


dr'^ 


1/2 


A\^ 

dr 


(33) 


where 


1 

P+c 


1 + e p (^ -L 1 + £ /i') 


A 


- A = o. 

PiG 


(40) 


Equation (HOj) corresponds to the cubic equation derived 
by Kagan [4/l| . 

Here we present a simple approximate solution to 
A. Since the bubble growth rate is significantly slower 
than the thermal velocity Vth, we assume that A <g; 
3aVth/i’c/(l + e)- Then, the factor in Eq. (l40l) is rewritten 


A = 


■ d 

( dr \ 



as 



U)J 

Tc 

(34) 




The growth rate of a spherical bubble is described by 
the Rayleigh-Plesset equation |27| |. 


3a Vth 


3aVth 


r.A 


(41) 


d^r 




dP 


m 


2'y Arjdr^., 

— 35 

r rdt 


which reduces Eq. (gni) to a quadratic equation and we 
obtain an approximate solution of A as 

A = Aine 


where 77 is the viscosity of the liquid. The evaporation 
rate from the bubble surface, i.e., the time evolution of 
the molecule number in a bubble i, is given by 0,111 



A- A- 

-^ine -^ine 

' 9 A • 

= 1 , ^-^VIS 


p — p 

JTeq -^g 


di , 9 a 

— = 47rr -- . , 

dt 1 + e y/27TmkT 


(36) 


+ 1 - 

A 

-^ine 

-^ine I 

d 2 ) 

94 

2 Avis I 

A . 

-^vis: 

and .Ajjje 

are deter- 


mined by the evaporation, the inertia of the fluid motion, 
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and the viscosity at the surface region of the bubble, re¬ 
spectively, and given by 


with 


Ae, 


A, 


27 aVth 

r2 (l + e)P3q’ 

7 

2rcTj’ 


27 \ 

pir^aj 


1/2 


(1 + £)^Peq 


(43) 

(44) 

(45) 


(46) 


In the limit Aeva ^vis, 4li„e, the growth is regulated 
by the evaporation process and A = ^evap- If ^ine (^vis) 
is the smallest of them, the growth is determined by 
the inertia (the viscosity) of the fluid and A is given by 
Aine (^vis)- We recall that this approximate solution of 
Eg. di^ is valid in the case where A <C 3at4th/?'c/(l + e)- 
To evaluate Jo of Eq. we also need (PAG/di'^ and 
di/dr. Using Eq. (15^ . we obtain 


d^AG r d'^AG \ r d^AG \ fdPA 
dr2 5^2 V 9PI )^\dr ) 


= —Sttj 



3Peq 

feq — -Pl 



Aeva. A 
Peq-Pl A ) 


(47) 


^ = 47rr2:^ / (1 + 
dr V 3Peq Aeva/ 


(48) 


Their Peq-dependencies affect Jq. Thus by the use of A, 
we find the expression for the prefactor in the nucleation 
rate: 



21 - 1/2 


(49) 


Note that the explicit expression of the exact prefactor 
was not presented by Kagan [4ll |. 

Fig. 2 shows the prefactor given by Eq. as a func¬ 
tion of the liquid pressure Pi for Lennard-Jones liquid 
with T* = 0.6,0.7, and 0.855. The equilibrium pressure 
and the surface energy are given by Peq = 0.0034etT“^ 
and 7 = 0.51ecr“^, respectively at T* = 0.6 [ 2 ^. Here, 
the evaporation coefficient a is set to unity and e to zero. 
The viscosity 77 is set to be 0.6 ertT“^ with the time unit 
T according to Angelil et al. [^ . 

When the liquid pressure is very close to the equilib¬ 
rium value, Jo agrees with the CNT value of Eq. (|4]). In 
this case, the bubble growth is regulated by the evapora¬ 
tion, and the mechanical equilibrium is almost satisfied 
on the path of bubble growth near the critical size. For 



FIG. 2. (Color online) The prefactor obtained by Eq. (l49l) 
as a function of the liquid pressure Pi. Panel a shows 
the case of T* = 0.6 (see text for detail). The solid line 
presents the prefactor with the exact solution of A to Eq. (PH) 
whereas the dashed line uses the approximate solution of 
Eq. (gH). Eg. (|42|l reproduced very well the exact solution 
within the accuracy of 0.1 %. The prefactors calculated with 
A — Aeva,Aine, or Avis are also plotted with gray dashed, 
light blue dotted, or green dotted-dashed curves. The liquid 
pressures used in the molecular dynamics simulations (see Sec 
III) are plotted by circles on the horizontal axis. The pref¬ 
actors for various temperatures T* — 0.6, 0.7 and 0.855 are 
shown in panel b. 


a large negative pressure, on the other hand, the prefac¬ 
tor can be significantly smaller than the CNT value by a 
factor ^ 10 ^ because of the viscosity and inertia effects 
in the liquid. In this case, chemical equilibrium holds 
rather than the mechanical one on the growing path. We 
obtain this growing path near the critical size in the r — i 
plane, by evaluating di/dr with Eq. (14811 . The obtained 
path is plotted in Fig. Id. 
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FIG. 3. (Color online) Comparison between the improved 
formula and the CNT for (TSF) Lennard-Jones fluid at T* = 
0.855. Panel a shows the comparison in J. The improved nu- 
cleation rate given by Eqs. (1281) and (1491) without the Tolman 
correction is plotted with the gray line. The improved formula 
with the Tolman correction is the thick red curve, where St 
is set to be 0.28a. The nucleation rate given by Eg. (1281) with 
Jo = To, CNT is the green dashed curve. Panel b shows the 
free energy for formation of a critical bubble for the PCNT 
without and with the Tolman correction (gray and red solid 
curves) and for the CNT (dotted curve). Panel c shows the 
critical bubble radii (same models and line-styles as above). 
In all panels, the results obtained by the molecular dynamics 
simulations are plotted with filled circles. 


III. COMPARISON BETWEEN THE 
CLASSICAL NUCLEATION RATE AND THE 
EXACT FORMULA 

We first compare the theoretical nucleation models de¬ 
scribed in Sec II with the molecular dynamics simulations 
by [ 2 ^ . The MD simulations used molecules with a trun¬ 
cated force-shifted (TSF) Lennard-Jones potential and a 
cutoff length of 2.5cr. In the models, we thus use the val¬ 
ues of TSF Lennard-Jones liquid for the thermodynam¬ 
ics data such as the surface tension and the equilibrium 
pressure. The evaporation coefficient a is assumed to 
be unity. In the calculations of the improved nucleation 
rate, we use Ea. d^ . In Figs. El we show comparisons 
at T* = 0.855. We plot the results of MD simulations 
(with filled circles) and four theoretical models. The first 
model is the CNT. The CNT nucleation rate is given by 
Eq. (El). Other nucleation rates are normalized with re¬ 
spect to this in Fig. Eli. The second (green dashed curve) 


is the PCNT nucleation rate which is given by Eq. (l28l) 
with Jq = Jo, CNT of Eq. ([1]). The Poynting correction 
factor is evaluated from Eqs. and (E5)) to include the 
non-ideal gas effect. 

The third (gray curve) is our improved expression given 
by Eqs. (1^ and ((1^ . The last one (red curve) uses the 
same improved expression as the third one, but also in¬ 
cludes the Tolman correction to the surface tension. The 
Tolman correction describes the size-dependent surface 
tension as 

7 = 7oo/(1 + 2(5t/t), (50) 

where 700 is the surface energy of the planar interface. 
The Tolman length in the Tolman correction St expresses 
the curvature dependence on the surface tension of the 
bubble. If St is positive, the surface tension of small 
bubbles or droplets is smaller than the planar one. 

In Fig. El^i we find that the improved formula (gray 
curve) can be smaller than the CNT (the PCNT) by a 
factor of ~ 10“^ (~ 1/5). Fig. El shows that the Tolman 
correction is necessary to reproduce the results of the 
MD simulations. From the fitting, we obtain St = 0.28a. 
Fig. Eb shows the peak values of the free energy for bub¬ 
ble formation. We find the PCNT correction almost re¬ 
produces AG(rc) in the MD simulations. Fig. Et shows 
the critical bubbles, Tc, obtained from the theoretical 
models and MD simulations. All models successfully re¬ 
produce the Tc measured in the MD simulations. Note 
that the Tolman length of 0.28cr obtained by fitting is 
much smaller the critical radius. Thus the correction to 
the surface tension of Eq. (1501) contributes only at the 
10 % level, yet significantly affects the nucleation rate 
due to the high sensitivity that the surface tension has 
on J (see Eq. (E51) ). 

Fig-Slshows the comparisons in J at T* = 0.6,0.7, and 
0 . 8 , where the equilibrium pressure and the surface en¬ 
ergy are given by Peq = O.OSOSecr”^ and 7 = 0.17ea~‘^, 
respectively at T* = 0.8 [ 2 ^. We find that deviations 
in the improved formula (gray curves) from the CNT 
can be several orders of magnitude in all cases. Com¬ 
pared with the values obtained by MD simulations, the 
improved formula without the Tolman correction gives 
significantly lower nucleation rates. By the fitting, we 
obtain St = 0.29cr, 0.32cr, and 0.38cr for T* = 0.8, 0.7, 
and 0.6, respectively. The fits at T* = 0.6 — 0.855 show 
that the Tolman length depends weakly on temperature. 
Since the radii of the critical bubbles are 3-7 a in the 
MD simulations, we see from Eq. dsni) that the obtained 
Tolman lengths of ~ 0.3fT correspond to corrections of 
10-20 % to the surface tension from the bulk values. As 
seen in Eigs.ElandlU these small corrections to the sur¬ 
face tension actually improve much the predictions of the 
nucleation-rate formula at T* = 0.6 — 0.855. 

We also make a comparison with the argon bubble nu¬ 
cleation experiments by [l^. Fig. E] shows the compar¬ 
isons in J at T = 137 K and 110 K. As for the thermody¬ 
namics data in models such as the surface tension and the 
equilibrium pressure, we use the data of argon [i^. We 







FIG. 4. (Color online) The same as Figl3Ta). but for T* — 0.8, 
0.7, and 0.6. In the improved formula with the Tolman cor¬ 
rection (red curves), the Tolman length 5t are set to be 
5t = 0.29, 0.32, and 0.38 a for T* = 0.8, 0.7, and 0.6, re¬ 
spectively. 






FIG. 5. (Color online) Comparison between the theoreti¬ 
cal models and laboratory experiments in J for argon at 
T = 137 K (a) and 110 K (b). The improved nucleation rates 
without the Tolman correction and with the Tolman correc¬ 
tion of 5t = 0.3cr at r = 137 K {5t = 0.13ct at T = 110 K) 
are shown by the gray and red solid curves, respectively. The 
nucleation rate given by Ea. (l28ll with Jo = Jo.evap is shown 
by the green dashed curve. The CNT is shown by the dot¬ 
ted curve. In the thermodynamics data such as the surface 
tension and the equilibrium pressure, we use the data of ar¬ 
gon [^. We also show the experimental results [T^l (circle), 
where the error in the liquid pressure is determined from the 
two points at T = 137 K in [l^ . 


find the deviations between the CNT and the improved 
formula can be several order of magnitudes in these cases, 
too. 

The experimental result is consistent with all of the 
theoretical predictions within the error at T = 137 K. 
Thus it is difficult to fix the Tolman length in this case. 
For T = 110 K, the improved formula gives significantly 
smaller nucleation rates than the experimental result. 
The fitting with the Tolman correction indicates that 
the Tolman length 5t is 0.13(7 ± 0.14(7. The error in 
5t comes from the errors in J and Pi in the laboratory 
experiment. This value of the Tolman length can be 
marginally consistent with the ones from the fits to 
the MD simulations. To further constrain 6t more 
experimental data are needed. 


In Fig. 6, we plot the Tolman length for bubbles ob¬ 
tained from our analysis (red circles) as a function of T* 
as well as those in the previous studies (other red sym¬ 
bols). As for our fitting data with the argon experiment 
in Fig. 5b (red open circle), the normalized temperature 
T* is obtained as 0.92, because e/k = 120 K for argon. 
Our results indicate that the Tolman length is almost 
constant with temperature. Our results agree quite well 
with the previous results by 371. 

Baidakov and Bobrov [2^ and Blockjs^ obtained 
small negative values in their MD and MC simulations. 


respectively. The deviation in the Tolman length between 
their results and ours comes from the different definitions 
of St (or the different curvature dependence of the surface 
tension). Their obtained surface tensions are always less 
than the bulk values, which is consistent with our cur¬ 
vature dependence with a positive Tolman length. For 
example, Baidakov and Bobrov [2^ gave the curvature 
dependence of the surface tension as 

7 = 7oo/(l + 2(5t/?' + (51) 

If Eq. (15(1)) is used for the evaluation of St instead of their 
curvature dependence, their results give positive St- In 
fact Eq. dsni) with St = 0.3(7 also very successfully re¬ 
produces the surface tension measurements by Baidakov 
and Bobrov [2^ (in their figure 9). This means that the 
results by Baidakov and Bobrov also agree well with ours. 

Moody and Attard also obtained negative Tolman 
lengths at T* > 1 from their MC simulations. Even at the 
high temperatures, nevertheless, their obtained surface 
tension increases with the bubble radius at r (5 2(7 (see 
their Eigure 10 and 11), which is also consistent with our 
curvature dependence. 

Recent calculations with density functional theory [sol . 
Isijl also result in negative Tolman lengths, by using a 
different curvature dependence similar to |25| |. However, 
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FIG. 6. (Color online) The Tolman length for bubbles as 
a function of temperature obtained by this study and the 
previous ones [^ . (red symbols). Our analysis 

with MD simulations (filled circles) and the experiment at 110 
K (open circle) suggest the Tolman length is approximately 
given by St = 0.3 a. For reference, the Tolman len gth for 
droplets are also shown with gray symbols [^. [3^ - [3^ . IdOU . 


they also obtained smaller surface tensions than the bulk 
and, in this sense, their results are consistent with ours. 

For reference, we plot the Tolman length for droplets in 
Fig.m We find the value for droplets by [i^ is consistent 
with our results for bubbles, although there is some scat¬ 
ter in the previous values. Our results suggest that the 
Tolman length is approximately given by ~ 0.3 ct, and is 
temperature independent. The Tolman correction with 
St = 0.3 <7 significantly improves the prediction of the 
nucleation rate with our nucleation rate model. 


processes (i.e., evaporation, thermal conduction, viscos¬ 
ity, and inertia of liquid). This difference in the prefac¬ 
tor between the bubble and droplet cases has been over¬ 
looked in many studies. We have also compared the im¬ 
proved expression of the nucleation rate with results of 
the MD simulations and laboratory experiments. Our 
findings are summarized below. 

• In bubble nucleation, the prefactor is strongly de¬ 
pendent on the degree of non-equilibrium, i.e., the 
liquid pressure, whereas it is constant in the droplet 
case. In the case of a large negative liquid pres¬ 
sure (or the highly viscous case), the prefactor can 
be far smaller than the droplet case by a factor of 
10-3 - 10-1 (Fig. 2). 

• When the liquid pressure is slightly below the equi¬ 
librium pressure, the deviations in the free energy 
for bubble formation from the CNT become large. 
For example, in dark matter detection experiments 
which use superheated liquids as targets , the 
degree of non-equilibrium is very small. For such 
a near-equilibrium case, the PONT should be used 
instead of CNT. 

• Comparisons of our improved expression for the nu¬ 
cleation rates J to results from MD simulations 
and laboratory experiments suggest that the sur¬ 
face tension depends on the bubble size at the nano¬ 
scale level. The improved expression including the 
Tolman correction to the surface tension with a 
small Tolman length of ~ 0.3 (t leads to good agree¬ 
ments with the recent MD simulations and labora¬ 
tory experiments (Fig. 6). 
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